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Abstract 


A study of delamination is performed 
including strength of materials and fracture 
mechanics approaches with emphasis placed on 
methods of delamination prediction. Evidence is 
presented which supports the inclusion of the in- 
plane stresses in addition to the interlaminar stress 
terms in delamination criteria. The delamination 
can be modeled as a resin rich region in between 
ply sets. The entire six component stress state in 
this resin layer is calculated through a finite 
element analysis and inputted into a new Modified 
Von Mises Delamination Criterion. This criterion 
builds onto previous criteria by including all six 
stress components. The MVMDC shows 
improved correlation to experimental data. 


Problem Scope 


The problem of free edge delamination in 
composite laminates has been studied for many 
years. While much progress has been made, there 
still exists the need for simple tools with which 
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designers can evaluate candidate laminates for 
their propensity to delaminate. Such a tool could 
improve applications where delamination failures 
are either expected or would be very detrimental. 
Delamination criteria have been proposed in 
recent years, however this investigation will 
explore an alternate approach. 

There are two main routes taken when 
studying the free-edge delamination problem: 
strength of materials and fracture mechanics. The 
strength of materials approach defines failure to 
occur when a critical stress state is reached based 
on material allowables. On the other hand, the 
fracture mechanics approach suggests that the 
initiation and growth of a crack depends on the 
energy state at that local station. In the present 
study delamination initiation is considered the 
failure of the part. The growth of the delamination, 
stable or unstable, is the continuation of that 
failure and is not addressed. 


Previous Work in Delamination Criteria 


Among the first to research the problem 
were Puppo and Evenson' who studied 
delamination by modeling anisotropic layers 
separated by thin isotropic shear layers. They 
recognized the stress conditions and laid the 
foundation for further study. Pipes and Pagano 1 2 
followed with a series of studies over a number of 
years. Their efforts dealt mainly with the 
mechanics of the free edge problem and solution 
techniques for the stresses including finite 
difference, elasticity, and global-local solutions. 
They also investigated the effects of laminate 
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thickness and stacking sequence. They 
demonstrated the existence of interlaminar 
stresses in a boundary layer near the free edge of 
a laminate. This boundary layer is suggested to be 
on the order of the laminate thickness. These 
interlaminar stresses can become large and lead 
to delamination. The interlaminar stresses and 
coordinate orientation are shown in Figure 1. 
They also concluded that the sign and magnitude 
of these stresses are functions of the material 
properties, fiber directions, thickness and stacking 
sequence. The stress condtion as calculated by 
Pipes and Pagano is shown in Figure 2. Pipes 
and Daniel* provided experimental verification of 
Pipes and Pagano results. They also concluded 
that a three-dimensional stress field arises in a 
boundary layer of approximately one laminate 
thickness from the free edge. 



Figure 1 : Laminate Coordinate System and Stress directions 
(from Jones”’) 

Isakson and Levy" performed some early 
finite element analyses by modeling plies as 
homogeneous, orthotropic membranes separated 
by isotropic shear layers which only develop 
interlaminar shear stresses. They calculated the 
interlaminar stresses and found decent agreement 
with existing values. Rybicki' 2 extended the 
procedure to a three-dimensional analysis which 
achieved better agreement with Pagano solutions. 
The in-plane stresses at the ply interface were 
also calculated at the edge and compared with 
values obtained from the middle of the laminate. 

As the stress state at the free edge 
became more characterized, some researchers 


began to suggest that the interlaminar stress state 
became singular at the edge. This was 
investigated by Wang and Crossman' 3 who studied 
the singular behavior through finite element 
methods with fine mesh resolution. They 
calculated significant in-plane stresses as well as 
interlaminar stresses and suggested that these 
singular stresses dissipated into the material 
through a stress relaxation that degrades the 
material in that local region. Hsu and Herakovich' 4 
provided mathematical evidence of singular 
interlaminar shear stresses and a finite 
interlaminar normal stress. Wang and Choi' 5 
performed an extensive analysis on the edge 
singularity and concluded that the state of stress is 
inherently three-dimensional and that the 
interlaminar shear term exhibits a weak singularity 
which is dependent only on the material elastic 
constants. 



Figure 2: Interlaminar stresses as computed by Pipes and 
Pagano (from Pagano*) 


As the development of the stress state at 
the free edge matured, the treatment of these 
possibly singular stresses in a quantitative sense 
became important. The need to address this 
mathematical behavior with a physical quantity 
was required. Whitney and Nuismer'* introduced 
the idea of an average stress in notched 
composite laminates. Their model averaged the 
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longitudinal stress over a certain distance and 
predicted failure when this average stress 
exceeded a critical value based on material 
properties. This concept of the average stress 
was applied to the free edge by Kim and Soni' 7 " 
who averaged stresses over an arbitrary critical 
length of one ply thickness from the free edge. 

Many strength of materials models for 
delamination described thus far in this paper have 
focused on the interlaminar shear and normal 
stresses as the sole contributors to delamination. 
Much debate is given to the extent which either or 
both the interlaminar normal stress, o n , and the 
interlaminar shear stress, a a , contribute to 
delamination. Also debated is how the sign of 
these terms changes their contribution. Early 
speculation by Pipes et al' 9 suggested that 
delamination is a complex condition where 
stresses and failures interact. This notion of 
interaction was submitted by Kim and Soni 1 ' with 
the interlaminar normal and shear stresses. While 
the in-plane stresses have been calculated as 
significant, they have not been considered as 
playing a role in delamination. Traditionally the in- 
plane and out-of-plane stresses have been placed 
into separate criteria. This is primarily due to the 
consideration of delamination as the failure of an 
interface between plies rather than a region 
between plies. 

In a similar development, fracture 
mechanics has been applied to the delamination of 
composite materials. The underlying premise is 
that a delamination is a crack in the thin, isotropic, 
homogeneous, and brittle matrix material between 
two layers of fibers. There is a distinction between 
delamination and transverse cracking. Several 
researches who pursue the application of fracture 
mechanics to composite materials include 
transverse cracking as well as delamination in 
their studies. Fish and Lee 20 demonstrated how the 
strain energy release rate significantly increases 
with the onset of transverse cracking. In addition, 
Reifsnider et al 2 ' suggested that delamination 
interacts with transverse cracking in 90 degree 
plies. The present study considers transverse 
cracking a different failure mode and an attempt 
was made to only investigate laminates which do 
not show a propensity for transverse cracking. 

Rybicki 22 was among the first to discuss 
how a delamination could grow by the release of 
strain energy in sufficient quantity to create new 
surface area. This strain energy release rate or 


fracture toughness was investigated through finite 
element methods. Wang and Crossman” and 
Crossman et al 2 ' studied the application of the 
strain energy release rate, G c , to delamination both 
analytically and experimentally. They assumed 
that the strain energy release rate for delamination 
and transverse cracking were the same. S.S. 
Wang 25 studied mixed-mode stress intensity 
factors and G c terms. His findings show that for 
angle ply laminates delamination is a three- 
dimensional effect. An extensive examination of 
the application of fracture mechanics to 
composites is by O'Brien” 27 . The fracture 
mechanics approach is justified by the nature of 
the singular stress state which precludes the 
comparison with any strength parameter. He 
developed a simple expression for the strain 
energy release rate by comparing the difference in 
strain energy between a delaminated portion and 
an undelaminated portion of a laminate by 
comparing the two modulii. This relation can be 
used to determine a critical value of strain energy 
release rate and can be used to predict 
delamination. Assumptions of delamination size 
and strain energy are required. O’Brien also 
concluded that interlaminar flaws at the free edge 
either initially in the material or caused by a stress 
singularity do not cause delamination until a critical 
strain energy release rate is reached. 

In recent years, Brewer and Lagace 2 * have 
developed a criterion under the strength of 
materials approach which utilizes average 
stresses compared with laminate strength terms. 
The stress state was calculated through the force 
balance method developed by Kassapoglou and 
Lagace 29 . Their Quadratic Delamination Criterion 
(QDC) predicts failure based on critical values of 
both interlaminar shear and normal stresses. The 
QDC is a two parameter model in which the 
averaging dimension and interlaminar shear 
strength serve as fitting parameters. In some 
cases the interlaminar shear was held constant 
and the averaging dimension was used solely as a 
fitting parameter. Their analytical work was 
correlated to a series of experiments by comparing 
predicted values of delamination onset from both 
the QDC and strain energy release rate method to 
measured data. It is the intent of this work to 
extend the concept of the QDC to include all stress 
components and provide comparable agreement 
with the data. 
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Free Edge Stresses 

At the free edge of a flat laminate loaded 
in axial tension there arises a six component 
stress state which contains significant interlaminar 
shear and normal terms. These stresses become 
large within a boundary layer which has been 
suggested is on the order of the laminate 
thickness". These interlaminar stresses arise 
from mismatches in elastic properties between the 
plies. The laminate coordinate system, stress 
designations and the stress distribution are shown 
in Figures 1 and 2. 

The calculation of these stresses has 
been investigated for many years as was 
described earlier. In an effort to develop a design 
tool it is desirable to simplify these calculations as 
much as possible. Modern computer methods 
have enabled this through simple and efficient 
finite element methods which were employed for 
the calculation of the stress state. 

When considering the full stress state at 
the free edge, it is important to recognize the 
continuity of each of the components. The 
interlaminar stresses are continuous through the 
thickness at the free edge, but the in-plane 
stresses are not. In a [0/90], laminate the 

difference in the in-plane stresses in the two plies 
is obvious. At the interface between two ply sets 
an additional layer is considered which represents 
the resin rich region. This region is an artifact of 
the manufacturing process and is the site of 
delamination investigation. Through finite element 
methods, it is possible to obtain the six component 
stress profile in the region and subsequently input 
them into a criterion for delamination prediction. 

A finite element analysis was employed 
which captured both the material and geometric 
considerations of the edge of a flat laminate in 
axial tension. The stress state at the free edge 
was determined using a quasi-three-dimensional 
finite element approach with a 12 degree-of- 
freedom assumed-stress hybrid element based on 
the Hellinger-Reissner principle 3031 . A constant 
strain loading was applied which is continuous 
through the thickness. The model consisted of 
elements both through the thickness and along the 
width of a laminate. A resin layer was modeled 
between each of the ply sets. This layer allowed 
for the direct calculation of the in-plane stresses 


which are not continuous through the thickness at 
the edge. The full six component stress state was 
calculated in the resin region. A closeup of the 
edge of a typical model is shown in Figure 3 for a 
simple four ply laminate, [15/- 15],. The total width 
of the model extends to a large number of ply 
thicknesses to ensure the recovery of classical 
laminated plate theory in the center of the 
laminate. The element size at the edge is one 
quarter ply thickness in width and increases in size 
away from the edge. The element size in 
thickness is constant throughout the model and is 
one half ply thickness. The resin layer element 
thickness is one tenth ply thickness. 



Figure 3: Finite element mesh at the edge of a [±15], 

laminate. 


Figure 4 shows the profile of all six stress 
components near the free edge. The stresses 
obtained show excellent agreement with those 
obtained by other researchers using other 
methods 8 . Some additional runs were made to 
investigate the effect of element size on the 
average stress levels. Clearly, the singularity at 
the edge yields an increasing interlaminar shear 
value with decreasing element size. In the finite 
element analysis, the element size was kept to a 
practical level. The arbitrary averaging dimension 
was selected to be on the order of the laminate 
thickness and the average stresses over this 
distance were examined for various element 
resolutions. It was observed that one quarter ply 
thickness elements in width sufficiently captured 
the stress levels while remaining a practical size. 
Other iterations with finer meshes did not produce 
a significant change in average stress values. 
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Figure 4: Stress distribution at the edge of a [±1 5], laminate. 


It was concluded that the effect of the a 2l term is 
negligible and that the compressive o„ term is 
either negligible or predicts an initiation value 
much greater than in-plane failure values. 
Therefore, only tensile values of interlaminar 
tension were considered. The following form was 
used for correlation to the experimental data: 



Several researchers have concluded that the 
interlaminar shear term dominates delamination 
while others postulate that the tensile interlaminar 
normal term gives rise to the failure. This criterion 
contains both terms, however delamination is still 
attributed solely to out-of-plane stresses. 

The averaging dimension was optimized 
to fit the data best at 0.178 mm or about 1.33 ply 
thicknesses. The interlaminar shear strength was 
assumed to be 105 MPa which is consistant with 
values for this material system found in the 
literature. This value, as suggested by Kim and 
Soni' 7 is the in-plane shear strength of the 
material. 


Quadratic Delamination Criterion ” 


Proposed Criterion 


As was mentioned earlier, the Quadratic 
Delamination Criterion indicates delamination by 
comparing interlaminar shear and normal stresses 
with appropriate laminate strengths. The stresses 
used are average values computed over an 
arbitrary averaging dimension. The QDC is of the 
form: 



The proposed Modified Von Mises 
Delamination Criterion (MVMDC) builds on the 
QDC by retaining all stress components and ply 
strength values and is of the following form: 


gf-g.Os+q | | + 

R'R c + R'R C 


■(M 


z' 


-c N* 
Or 


where: 


where: 

Z*' = Interlaminar shear strength 
longitudinal shear directions 
Z* 7 = Interlaminar shear strength in 
transverse shear directions 
T = Tensile interlaminar normal 
strength 

T = Compressive interlaminar normal 
strength 


a|,an= In-plane principle stresses in the 
resin layer 

a M = Vectorial sum of the interlaminar 

shear stresses, o„ and 

0 ^ = Average interlaminar normal stress 

R 1 = Tensile strength of resin 
R c = Compressive strength of resin 
Z* = Interlaminar shear strength of ply 
Z* = Interlaminar tensile strength 
r = Interlaminar compressive strength 
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The stress terms are shown in Figure 5 
and the following paragraphs describe each of the 
terms in the above criterion: 

The first term contains the in-plane 
stresses. At the free edge, the in-plane stresses 
are not zero and contribute to the stress condition 
and should be included. Their inclusion was 
treated from a von Mises perspective. The 
average in-plane stress terms o„, and a ^ are 
rotated into principle stresses o, and o„. These 
stresses can be considered in a von Mises 
relation; however the characteristic strength c 0 
may not be a traditional strength parameter for 
anisotropic materials. In addition, the resin layer is 
an epoxy system in this investigation which may 
exhibit different strengths in tension and 
compression. An examination was performed by 
Williams” where the von Mises criteria was 
applied to polymers exhibiting this effect. This 
resulted in the separation of the characteristic 
strength into the product of tensile and 
compressive strengths. In a laminate, these are 
assumed to be the transverse tensile and 
compressive strengths, Y and Y*. The in-plane 
term takes the form shown above. The three- 
dimensional von Mises failure theory for isotropic 
materials was applied to the stress state in the 
resin, however it did not give satisfactory results. 
This can be attributed to the layer effect because 
of the thin layer in the z direction. 

The second term is the interlaminar shear 
term. It is the vectorial sum of the average 
interlaminar shear stresses and In many 
cases one of the two terms will significantly 
dominate the result, however they are both 
included. Assuming transverse isotropy the same 
strength value for interlaminar shear can be 
compared to both terms. Alternately the stresses 
can be considered individually as in the QDC. 

The third term is the interlaminar normal 
term. It accounts for both tensile and compressive 
interlaminar normal stress contributions to 
delamination. It has been suggested in the 
literature that compressive interlaminar normal 
stresses do not contribute to delamination. In 
some cases these stresses were thought to inhibit 
delamination, but there is no physical evidence to 
suggest this. The average stress is compared to 
the interlaminar normal strengths. It is important 
to note that the a u term often exhibits a crossover 
point from tension to compression within a couple 


of ply thicknesses. Care should be taken to 
consider this when computing an average stress 
as it is possible for the stress state to average out 
to zero when actual tensile and compressive 
stresses exist within the region. This effect was not 
witnessed in this study. While the interlaminar 
tensile strength can be obtained though flatwise 
tension testing, the compressive value is difficult to 
obtain. Current experimental methods may not 
accurately predict this value requiring it to be 
estimated. These values are assumed as the 
transverse tensile strength, Y, and the transverse 
compressive strength, Y". 





Figure 5: Stress directions for MVMDC 


The QDC excludes terms because of their 
negligible impact on the resulting prediction, 
however there has been no physical evidence 
presented that some components should be 
excluded. It is often difficult to intuit how 
significant a stress value will be prior to the 
detailed stress analysis. In addition, modern 
analytical methods are frequently coded in 
computer programs which easily perform the 
additional computations. Therefore the entire six 
component stress state is represented in the 
MVMDC. 

The validity of the MVMDC is shown by 
both applying it in a uniaxial configuration and by 
comparing it to in-plane criteria. An additional 
finite element analysis was run with all zero 
degree plies and the stresses were backed out for 
the thin resin layer in between two adjacent plies. 
When the MVMDC is applied to these results, a 
strength is predicted which is larger that the in- 
plane strength of the laminate as predicted by in- 
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plane criteria and experiment. The MVMDC does 
not indicate any failure prior to in-plane lamina 
failure. 


Figure 6 shows the failure envelope for 
the MVMDC as compared to the Maximum Stress 
and Tsai-Wu criteria 33 . The MVMDC envelopes 
the in-plane criteria except for the small region in 
the first quadrant. For a condition where Classical 
Laminated Plate Theory holds, namely an infinite 
width laminate and strain continuity through the 
laminate, the MVMDC reduces to a plane stress 
von Mises relation based on resin values and does 
not indicate in-plane fiber behavior. 




Figure 7: First quadrant of MVMDC failure diagram 


Correlation to Experimental Data 


The MVMDC was applied to the same 
data as the QDC and indicated delamination very 
well. Figures 8, 9, and 10 show the correlation of 
the MVMDC to the data presented with the QDC. 
The plots show the comparison of the predicted 
divided by the experimental value for each of the 
members of each laminate family. An accurate 
prediction has a value of unity while a number 
greater than one represents an overprediction. 
The laminate families are [±15] n , [±15/0] n , and 
[0/±15] n where n is one to five. The material 
system is AS1/3501-6 graphite epoxy with the 
following assumed strength parameters: 


Figure 6: Modified von Mises Delamination Criterion Failure 
Diagram. 


R' = t = Y = 54 MPa 
R c = T = Y* = 236 MPa 
r= 105 MPa 


The three plots are very close in the first 
quadrant, however the region is very small and is 
a artifact of the material parameters used. This 
effect changes with different materials and the 
small overlap may or may not show up with 
different material systems. The first quadrant is 
shown in Figure 7. Regardless, the MVMDC was 
not intended for in-plane failure analysis. Lamina 
criteria such as Tsai-Wu or Maximum Stress 
should be used to predict in-plane behavior. 


The values indicated by the MVMDC are 
in close agreement with the data. The averaging 
dimension was 0.402 mm or 1 .75 ply thicknesses. 
This is seven elements of one quarter ply 
thickness width. This is about a 30% increase in 
the averaging dimension from the QDC 
investigation, however it is still on the order of 
what researchers suggest 14 ' 7 . The averaging 
dimension was determined by selecting whole 
numbers of elements and comparing the resulting 
predictions to the data. This parameter is arbitrary; 
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there is no evidence suggesting that this value 
represents any physical constant, however it has 
been suggested it may be a material parameter' 7 . 
The averaging dimension is used in this analysis 
solely as a fitting parameter and was optimized 
only to the mesh resolution of the model. The 
strength values were held constant reducing this to 
a one parameter model. 



Figure 8: Data points for [±1 5]„ laminate family 



Figure 9: Data points for [±1 5/0], laminate family 



Figure 1 0: Data points for [0/±1 5], laminate family 


The MVMDC predicts delaminations for 
the [±15] n and [0/±15] n laminate families that is 
closer to the experimental values than those 
predicted by the QDC. The QDC predicts slightly 
closer values for the [±1 5/0]„ laminate family than 
the MVMDC, but the difference is very small. 

Correlation with other data in the literature 
shows excellent agreement with the MVMDC. An 
attempt was made to distinguish between 
delamination and transverse cracking. Some 
studies of delamination include specimens which 
show a propensity for transverse cracking 2434 (90 
degree plies). Kim and Soni performed 
experiments and obtained some data which did 
not utilize laminates with angle plies greater than 
30 degrees. The material system used is 
T300/934 with published values for ply properties 
and strengths* 338 . The averaging dimension was 
changed to 1 .25t for the best fit while the strength 
parameters were kept constant. The correlation to 
this limited data set is shown in Table 1 : 



Table 1 : Correlation to data in the literature 
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Conclusions and Summary 


A brief history of the study of delamination 
and interlaminar stresses is presented as a 
background for the development of a delamination 
criterion. The following conclusions were drawn. 

1. Delamination failure can be characterized as 
the failure of a resin rich layer or interphase in- 
between ply sets. Finite element methods can 
be used to model this resin layer and compute 
the full stress state. 

2. When considering the stress state in a resin 
rich layer in between ply sets, A full six 
component stress state is as valid as 
individual interlaminar stresses when 
considering the conditions for delamination. 

3. The Modified von Mises Delamination 
Criterion has improved delamination prediction 
over the Quadratic Delamination Criterion and 
strain energy release rate methods. Strength 
values used were nominal material properties 
and could be obtained through normal 
experimental methods. The single parameter 
is the averaging dimension which is a curve 
fitting term which can be arbitrarily chosen and 
experimentally fit. 

The proposed criterion can be easily 
applied when a full stress solution is generated. 
Finite element methods can easily provide the 
stress state for a variety of laminate 
configurations. This method can then be applied 
to evaluate candidate laminates for their 
propensity to delaminate and aid in the more 
efficient, durable and cost effective composite 
designs. 
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